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ABSTRACT 

We study the line profiles of a range of transition region (TR) emission lines observed in typical quiet Sun 
regions. In magnetic network regions, the Si IV 1402 A, C IV 1548 A, N V 1238 A, O VI 1031 A, and 
Ne VIII 770 A spectral lines show significant asymmetry in the blue wing of the emission line profiles. We 
interpret these high- velocity upflows in the lower and up per TR as the quiet Sun equivalent of the recently 
discovered upflows in the low corona above plage regions (Hara et aD2008l) . The latter have been shown to be 
directly associated with high- velocity chromospheric spicules that are (partially) heated to coronal temperature s 
and play a significant role in supplying the active region corona with hot plasma (iDe Pontieu et al.l 12009). 
We show that a similar process likely dominates the quiet Sun network. We provide a new interpretation of 
the observed quiet Sun TR emission in terms of the relentless mass transport between the chromosphere and 
corona - a mixture of emission from dynamic episodic heating and mass injection into the corona as well as 
that from the previously filled, slowly cooling, coronal plasma. Analysis of the observed upflow component 
shows that it carries enough hot plasma to play a significant role in the energy and mass balance of the quiet 
corona. We determine the temperature dependence of the upflow velocities to constrain the acceleration and 
heating mechanism that drives these upflows. We also show that the temporal characteristics of these upflows 
suggest an episodic driver that sometimes leads to quasi-periodic signals. We suggest that at least some of the 
quasi-periodicities observed with coronal imagers and spectrographs that have previously been interpreted as 
propagating magnetoacoustic waves, may instead be caused by these upflows. 
Subject headings: Sun: magnetic fields - Sun: chromosphere - Sun: transition region - Sun: corona 



1. INTRODUCTION 

One of the most overlooked but important problems in so- 
lar and stellar physics is the nature of the transport mecha- 
nism that is required to supply enough mass from the cool 
chromosphere to the hot corona to counter the effects of 
the slow, but steady cooling and draining of the corona o n 
timescales of order half an hour (Schrij ver & Zwaanl 12000). 
Ever since their discovery in the late 1800s when "verti- 
cal flame s" where first seen at the limb during total solar 
eclipse fll eccfii 1 87 7b, spicules or dynamic jet-like extrusions 
(Roberts 1945; Beckers fl968|) have received much attention 
as agents in the trans port of energy and mass in the lower so- 
lar atmosphere (e.g.,lThomaslll948l: IPneuman & Kopp|[T978l: 
Athav & Holzer 1982). This is because the mass they carry 
upward is estimated to be two orders of magnitude larger than 
that needed to supply m ass to the corona and solar wind (e.g., 
Athay & Holzer 1982). One shortcoming was that a signature 
of these ejecta could not be obs erved at "typical" coronal tem- 
peratures fe.g JWith broe 1983) and so the concept of spicules 
as a hot mass transport mechanism for the outer solar atmo- 
sphere has been largely mothballed. 

The Sola r Optical Telescop e (SOT; iTsuneta et aO 120081) 
of Hinode ( Kosug fet al.l l2007|) has revolutionized our view 
of the dynamic chromosphere r evealing the existence o f 
at least two types of spicules (|De Pontieu et al.l l2007b|) . 
"Type-I" spicules are long-lived (3-5 minutes) and ex- 
hibit longitudinal motions of the order of 20km/s that 
are driven by shocks resulting from the leakage of p- 
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modes and convective motions, as well as magnetic en- 
ergy release in regions around stro n g magnetic flux con- 
centrations ([b e Pontie uTt al.l 120041 LHa nsteen et al. 2006; 
IDe Pontieu eta L 2007a: iMartinez-Svkora et al.ll2009l) . Their 
contribution to TR emission may be limited to a hot sheath 
with lit tle of the chromospheric plasma heated to TR temper- 
atures (Ide Wijn & De Pontieu! 12005 1 Judgel 12008b . "Type-H" 
spicules, on the other hand, inhabit a truly differen t dynamic 
regime than the "classical" spicules studied by Roberts (1945) 
(or others) and have a direct impact on the TR. These recently 
discovered spicules show much shorter lifetimes (50- 100s) 
and higher vel ocities (50-150 km/s), a nd are likely caused by 
reconnection (|De Pontieu et al. 2007b). Type II spicules often 
exhibit only upward motion, with the feature fading along its 
entire length in a matter of seconds, which is suggestive of 
r apid heating to (at least) TR temperatures. 

be Pontieu eTall (120091) (hereafter D2009) recently studied 
active region plage and discovered a direct connection be- 
tween Type-II spicules and a persistent, faint, but strongly 
blue-shifted component of emission in hot coronal lines at 
the magnetic footpoints of loops. Their work suggests that 
coronal heating often occurs at chromospheric heights in as- 
sociation with jets that are driven from below and supply the 
corona with a significant fraction of the required mass flux. 

Here we show that a similar process occurs in the quiet Sun 
for a wide range of temperatures that cover both the lower and 
upper TR (or low corona). In contrast to D2009, the current 
paper is focused on the quiet Sun. We study several aspects of 
this mass-loading process that were not addressed in D2009. 
In § 2 we discuss the quiet Sun observations and present a 
novel ana lysis technique for U V emission lines observed b y 
SUMER (Wilhelm et al.lfl995h on SOHO dFleck et al.lll995h . 
We perform a detailed study of the strengths and limitations 
of the line asymmetry analysis. We investigate the effects of 
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blends on the analysis, and compare data observed in both 
first and second order to show that our conclusions are not 
impacted by blends. We describe the temporal behavior of the 
blueward asymmetries and compare them to those of type II 
spicules. Most importantly, we put forth a new interpretation 
of the nature of the emission from the quiet Sun transition 
region in terms of the mass cycling associated with the heating 
and cooling of the plasma associated with the type-II spicules. 
We discuss this new interpration in § 3. 

2. OBSERVATIONS & ANALYSIS 

In the following we will present line asymmetry analy- 
sis of SUMER raster scans for spectral lines of the follow- 
ing ions: Si IV 1402 A, C IV 1548 A, N V 1238 A, 
O VI 1031 A, and Ne VIII 770 A, which are formed over a 
wide range of temperatures (under equilibriu m conditions), 
respe ctively for logT = 4.8, 5.0, 5.3, 5.5, 5.8 dMazzotta et al.1 
1998). In addition, we show a sit-and-stare temporal sequence 
in N V 1238 A and analyze the episodic nature of the blue 
asymmetries. 

2. 1 . Simultaneous C IV 1548 A and Ne VIII 770 A 
observations 

The main data set used in the analysis presented is 
an absolutely calibrated SUMER spectroheliogram in the 
1530- 1555 A spectral range on September 22 1996 from 
00:40-08 :28UT for whi ch further details can be found in 
Hassl er et al.l (Il999|) and lDavey et al.l (|2006|) . This portion of 
the solar UV spectrum contains emission lines of singly ion- 
ized Silicon (Si II) with a rest wavelength of 1533. 43 A, seven 
times ionized Neon (Ne VIII) at 1540.84A (in second spec- 
tral order) and two lines of three times ionized Carbon (C IV) 
at 1548.20A and 1550.77A. The two C IV lines show iden- 
tical spectral behavior and so we only consider information 
from the former as it has better signal-to-noise. For reasons 
of brevity, low signal-to-noise, and the impact of line blends 
close to line center, we also exclude the Si II line from the 
analysis. This leaves us with two emission lines that span the 
TR, C IV and Ne VIII formed around 10 5 and 6xl0 5 K respec- 
tively. 

For context, Figure \T\ provides an EIT 
dDelaboudini ere et al.l l 19951) 195 A image (taken at 01:05UT) 
and a composite MDI (IScherrer et all [1995) line-of-sight 
(LOS) magnetogram 4 as well as the line intensities, non- 
thermal line widths (v n u determin ed from the recipe in 
iMcIntosh. De Pontieu & Tarbelll l2008h and absolutely cal- 
ibrated Doppler velocities of the C IV 1548 A and Ne VIII 
770.4A lines in the left and right columns respectively. The 
visibility of the supergranular network differs markedly 
between C IV (all network contours in Fig. Q] are based on 
C IV intensity) and Ne VIII (panel D) in which the network 
pattern is not obvious. We also see a correspondence between 
increased values of v nt and the LOS field strength 5 . In 

4 The magnetogram is a composite of the four 96-minute cadence full- 
disk magnetograms taken over the duration of the SUMER spectroheliogram 
(01:39; 03:15; 04:48; 06:27UT). The composite is needed to partially take 
into account the considerable evolution of the photospheric magnetic field 
during the 7 hours of SUMER observations. 

5 We note that there are many locations of enhanced v n t that straddle 
regions where th e magnetic field is expected to be significantly inclined, a 
topic covered by Mcintosh, Leamon & De Pontieu ( 2009), or the locations 
of "explosive events", see, e.g., Innes et al. (1997) and § 2.1. 
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FiG. 1. — Contextual information for this Paper. Top row: the EIT 195 A 
and composite MDI LOS magnetogram for the SUMER spectroheliogram of 
September 22 1996. Subsequent rows show the fitted line intensities, non- 
thermal line widths and absolutely calibrated Doppler velocities of the C IV 
1548 A (left) and Ne VIII 770.4 A (right) emission lines. In each case the 
fitted spectral information (panels C-H) are overlaid with the 150 count in- 
tensity contour of the C IV intensity as a proxy for the supergranular network 
boundary. 
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Fig. 2. — Reference network line profiles in C IV 1548 A (left) and 
Ne VIII 770.4A (right) for the network element centered on x=-373 // , 
y=452 // . Shown inset in each panel are the intensities, Doppler velocity, pro- 
file asymmetry at 50km/s (±12km/s; see text), v n t and MDI LOS magnetic 
field strength (in panel A only). The normalized mean profiles are shown for 
the network (blue solid line), the surrounding inter-network region (red solid 
line) while the Gaussian fits to those profiles and residuals (profile - fit) are 
shown as dot-dashed lines and triangular symbols of the appropriate color. 
The horizontal dashed line is the zero-level of the profile residual [placed at 
a normalized intensity of 0.2]. 

the bottom row of Fig. [T] the striking prevalence of red- 
shifted material is evident in the C IV emission (see, e.g., 
Gebbieetal. 1981, and many others) which is enhanced in 
the supergranular network vertices. In contrast, the Ne VIII 
emission in those same locations show modest blue-shifts 
or outflows (see Fig . [2]), that are typical for quiet Sun 
(IMcIntosh et al.ll2007h . 

We employ the technique discussed in D2009 to quantify 
the asymmetry of the SUMER line profiles. We first perform 
a single Gaussian fit (with a constant backgrou nd) to the line 
profile (IMcIntosh. De Pontieu & Tarbelll [200 8) to determine 
its centroid position. Then we interpolate the line profile to a 
wavelength scale that is ten times finer and calculate the to- 
tal number of counts in the interpolated profile in two narrow 
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range of wavelengths/velocities (two original SUMER spec- 
tral pixels wide, or 24 km/s at 1548 A) that are symmetri- 
cally positioned from the line centroid position, one in the 
blue and one in the red wing. The difference of the counts in 
the two wing positions, red-blue (R-B), provides a measure 
of the asymmetry of the profile at the chosen wavelength or 
velocity offset. Of course, a measure of zero indicates that 
the profile is symmetric in that velocity range. The R-B mea- 
sure is related to the third (skewness) and fourth (kurtosis) 
moments of the the line profile. For relatively noisy line pro- 
files, such as those studied, the difference of the two wings is 
a measure of asymmetry, provided that no spectral blends are 
present in the lines. 

We note that the Ne VIII 770.4A line, seen in second spec- 
tral order by SUMER, has t wo blends of S i I. As determined 
from the HRTS spectral atlas (Brekke 1993) these blends have 
equal magnitude, are ~20 times weaker, and are symmetri- 
cally distributed at ~20km/s either side of the Ne VIII rest 
position in the quiet sun. Based on the symmetry of the blend 
positions, their equal intensity and considerably lower ampli- 
tude (than Ne VIII 770.4A) we deduce that their net effect on 
the Ne VIII line R-B diagnostic in the quiet sun is negligible 6 . 
In a coronal hole the relationship between the three spectral 
lines is more complex and is discussed at length in a forth- 
coming Paper ( Mci ntosh. Leamon & De Po ntieu 2 009|) . 

Figure [2] illustrates normalized line profiles for a magnetic 
network element (red solid line) and the surrounding quiet re- 
gion (blue solid line) centered on coordinates (-373", 452"). 
We see that for both of the emission lines in the network ele- 
ment there is an excess of emission in the blue wing (similar 
to what is shown in D2009). The peak of the excess emission 
is at ~50km/s for C IV and ~55km/s for Ne VIII. The sur- 
rounding plasma shows little, if any, excess emission. This is 
clear from the plot of residual emission (triangular symbols 
that are raised by 0.2 from the lower portion of the plot) be- 
tween the mean region profile and a single Gaussian profile fit 
(appropriately colored dot-dashed line). In the network, the 
amplitude of the asymmetry is of the order of 3-8% of the 
peak brightness of the C IV line and 3-5% of the Ne VIII. 
Note that the excess emission in the blue wing occurs in a 
region in which the dominant core of the C IV spectral line 
predominantly shows modest redshifts of order a few km/s. 

Calculating the R-B asymmetry measure for the two emis- 
sion lines over a range of reference velocities from the line 
centroid, we produce the series of images shown in Fig. [3] 
(and accompanying movies). In both lines at low velocities 
(below sound speed expected for the relevant formation tem- 
peratures) we see that the R-B asymmetry is dominated by 
the red wing. This is particularly the case in the supergranular 
boundaries with the asymmetry becoming increasingly large 
over the strongest supergranular vertex flux concentrations. 
At higher reference velocities we see that the R-B asymmetry 
in the supergranular network switches sign indicating that the 
blue side of the line profile is becoming dominant. The change 
in sign is abrupt for C IV, and occurs at about 10-20km/s 
higher velocities for Ne VIII: we see a striking similarity be- 

6 In general, the presence of a strong blend close to the line core can impact 
the R-B asymmetry analysis because of its effect on line centroiding. For 
example, a strong blend in the red wing, close to the core of the line, will 
shift the line centroid slightly to the red, which can lead to a generalized 
weak blue asymmetry. Since the Si I blends are symmetrical, this is unlikely 
to be the dominant cause of the blueward asymmetries we find in Ne VIII 
770 A as observed in second order. 
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Fig. 3. — Red-blue (R-B) profile asymmetries (in arbitrary units of inten- 
sity) for the C IV 1548A (left) and Ne VIII 770.4A (right) emission lines. 
Again, for reference, the panels of the figure are overlaid with the 150 count 
intensity contour of the C IV intensity as a proxy for the supergranular net- 
work boundary. The electronic edition of the journal has movies showing 
the panels of this figure along with those at intermediate velocities and the 
components of Fig. [Q 
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Fig. 4. — Red-blue (R-B) profile asymmetries (in arbitrary units of inten- 
sity) corrected for the line centroid positions using the data presented in pan- 
els F & G of Fig. □respectively for the C IV 1548A (left) and Ne VIII 
770.4 A (right) emission lines. For reference the panels of the figure are over- 
laid with the 150 count intensity contour of the C IV intensity as a proxy for 
the supergranular network boundary. 
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Fig. 5. — Absolute velocities of the blue- wing asymmetry and related his- 
tograms in the C IV 1548 A (left) and Ne VIII 770.4 A (right) emission lines. 
The fractional intensities are computed for both emission lines everywhere 
the C IV intensity is greater than 50 (solid white contours). 




Fig. 6. — Fractional intensities and related histograms for the blue-wing 
asymmetry in the C IV 1548 A (left) and Ne VIII 770.4 A (right) emission 
lines. The fractional intensities are computed for both emission lines every- 
where the C IV intensity is greater than 50 (solid white contours). 

tween the two panels when we compare the R-B asymmetry 
for C IV at 50km/s and Ne VIII at 70km/s. 

We should note that the preponderance of red signal in the 
high velocity (> 70 km/s) panels of the left (C IV) column 
of the figure is the result of the analysis sampling the Si I 
1548.72A line that is 0.5A (i.e., 100 km/s) to the red of the 
C IV line center. This line is 4-5 tim es less bright than the 
C IV line (see, e.g. JCurdt et aLl feOOl) and we see that strong 
blue asymmetries still persist in many network locations, ie. 
indicating that the blue wing of the C IV line is still stronger 
than the narrow Si I line at the same velocity in the red wing. 
This blend reduces the visibility of the C IV asymmetry to 
higher velocities and can explain some, if not all of the veloc- 
ity offset between the C IV and Ne VIII asymmetries. We also 
note that a very weak Si I line occurs at about 100 km/s to the 
red of line center of Ne VIII 770 A (in second order). This 
line has an int ensity of order only 2% of the Ne VIII intensity 
(Brekke 1993). It likely plays a role in reducing the visibility 
of Ne VIII asymmetry at very high velocities (but since it is 
weaker than the blend in C IV it has much less effect). 

In order to place the derived R-B asymmetries in some rel- 
ative context we estimate their maximum amplitude and lo- 
cation in absolute wavelength space. This step requires that 
we combine the information presented in Figs. [T] and [3 us- 
ing the absolutely calibrated Doppler velocity from the line 
centroiding (panels QJj &QJI) to re-map the R-B asymmetry 
profiles. At each spatial pixel we adjust the velocity scale of 
the R-B profile to reflect the measured line centroid velocity; 
for example, a positive (red) line centroid Doppler shift will 
shift the entire R-B intensity profile to the red by the corre- 
sponding velocity and conversely for a blue Doppler shift of 
the line centroid. After shifting the R-B profiles we resample 
the panels of Fig. [5] to show the magnitude of the asymme- 
try in absolute velocity space, these are presented in the pan- 
els of Fig. |H In performing this analysis we can also extract 
the maximum fractional (relative to the peak line intensity) 
magnitude of the line asymmetry and its location in absolute 



velocity space. The plots of these quantities, and their corre- 
sponding histograms, can be found in Fig. [5] and Fig.0respec- 
tively where we see that the maximum fractional intensity of 
the asymmetry in the C IV and Ne VIII lines are of the order 
of 3-4% (increasing into the cell interior as the intensity falls 
off significantly) with velocities that are offset by ~20km/s. 

2.2. Ne VIII 770A first order observations 

To further illustrate the ubiquitous presence of these blue- 
wing profile asymmetries in the quiet Sun magnetic network 
we present data that does not suffer significantly from blend- 
ing with other lines. We also discuss the potential impact of 
low signal-to-noise on the detection of the subtle asymmetries 
we have found. 

We consider quiet Sun spectral atlas observations of the 
Ne VIII 770A line in first spectral order performed by 
SUMER close to the start of the mission (April 20 1997 
00:34UT: ICurdt et al.ll200Tl Fig. Ui and in a more r ecent ob- 
servation (July 3 2008 12:llUT7 lTia~n et al.1 [20091 Fig. 0). 
These data are taken with extremely long exposures (300s) of 
SUMER detector B and provide information about the Ne VIII 
770A line in first order. This line may be susceptible to Si I 
blends in the second spectral order, but is "clean" in the first 
order. We note that the nearest lines in the SUMER spectral 
atlas are the 769. 38A Mg VIII and 771.90A N III lines which 
have ~2% of the brightness of Ne VIII 770A line in the quiet 
Sun, a nd are at 390 and 590 km/s to the blue and red respec- 
tively (ICurdt et al.ll2001l) . 

Comparing Figs. [7] and [8] we are immediately drawn to 
the significant change in the strength, and "quality" of the 
SUMER spectra in the eleven years between the two observa- 
tions presented. The ratio of network-internetwork intensity 
(N/IN) across the entire spectral range is reduced by a factor 
of two, the N/IN profile is very noisy and of low spectral con- 
trast in the later observation, and the peak intensity of many 
of the emission lines (including Ne VIII 770A) in this spectral 
range are decreased by up to a factor of four. Indeed, com- 
paring the cross-slit variation of the Ne VIII 770A spectral 
quality measure (I / yT; a good estimate of the signa l-to-noise 
ratio - S/N, see, e.g. Chae, Schiihle & Lemaire 1998), we see 
a mean value of 11.8, with very little variance (1.5) and very 
little contrast between the network and inter-network regions 
in the July 2008 spectra. Compare these values to the April 
1997 case where the mean is ~15, with a variance of 4 and 
regions of the slit with values in excess of 20. Because the 
strongly blue- shifted component we are interested in is very 
faint (~ 5% of the line core), it is clear that the spectral qual- 
ity will significantly affect our ability to measure/quantify this 
component. Our analysis suggests that in the later phases of 
the SOHO mission SUMER spectra with a S/N of the order of 
20 (needed to make a robust detection of these faint emission 
components in the blue wing) are rare even in long exposure 
observations. This can explain the lack of blueward asymme- 
tries in the first order Ne VIII 77 A spectral profi les (using 
SUMER) as analyzed recently bv lTianet all (120091) . 

This issue of data quality is illustrated by the results of an 
R-B analysis of the SUMER Ne VIII 770A spectra in the ref- 
erence spectra of Figs. [7] and [8] A complicating factor in this 
region of the UV spectrum is the linear behavior of the con- 
tinuum background: even a small gradient can add additional 
weight to the red wing of the line unless it is taken into con- 
sideration when fitting the line center and when assessing the 
asymmetry of the profile. Therefore, we carefully fit a linear 
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Fig. 7. — The SUMER quiet Sun spectral atlas observation of April 20 
1997 00:34UT. The main panel of the figure shows several key measures: 
the SUMER reduced detector image displaying the spatial variation of sev- 
eral strong UV emission lines, the slit averaged spectrum (green), a vertical 
box enclosing the N IV 765A line. Isolating this line, we can compute the 
averaged TR brightness across the SUMER slit and identify sample regions 
of (supergranular) network and inter-network (horizontal boxes) labeled "N" 
and "IN" respectively. The red curve shows the Ne VIII 770 A I/Vl mea- 
sure as an estimate of the spectral quality across the SUMER slit. The white 
trace along the top of the detector image shows the ratio of the network to 
inter-network spatially averaged spectra over the spectral domain. The far 
right panel of the image shows the R-B analysis of the first order Ne VIII 
770A emission line along the SUMER slit. 
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Fig. 8. — The SUMER quiet Sun spectral atlas observation of July 3 2008 
12: HUT. The measures derived from the detector image are as for Fig. [7] 

background and perform the R-B analysis on the background- 
subtracted spectra. The results are shown in the far right pan- 
els of Figs. [7] and In the April 1997 reference spectrum we 
see a repetition of the signature observed in Figs.[5j[4l[5]and[6] 
(which were for the same spectral line, but in second spectral 
order): a red wing excess in the network regions at low veloc- 
ities transitioning to a blue wing excess at higher velocities. 
The gradual return to red wing excess at very high velocities 
is likely the result of the linear background which, when not 
fully removed far from the line core, will lead to a preference 
for red asymmetries. The fact that we observe blue asym- 
metries for velocities between 50 and 90 km/s despite this 
background issue underscores the strength of the blue wing 
excess 7 . In contrast, the R-B signature in the July 2008 spec- 
tra is highly dubious: there appears to be little correlation be- 
tween the patches of signal and any physical structure along 
the slit. We deduce that this erratic signal is the result of the 
considerably higher noise level of the spectra taken later in the 
mission. This degradation of the spectra is a subject that mer- 
its further investigation. Reference spectra like those shown 
have been taken periodically throughout the SOHO mission, 
so such an analysis is possible, but beyond the scope of this 

7 We note that the coronal hole reference spectrum of Ne VIII 770 A in 
first order (presented as the blue curve in Fig. 4, pg. 602, of Curdt et al. 2001) 
likely shows the blue-wing asymmetry as well. We suspect that the "knees" 
visible in the blue wings of the network/internetwork ratio in bright, spec- 
trally isolated, TR emission lines in the spectral atlas data (the green curve in 
Curdt et al. 2001) may be the signature of the network blue profile asymme- 
try. Further analysis and modeling are required to ascertain the validity and 
uniqueness of this point. 



Paper. 

2.3. Raster scans of Si IV 1402k, N V 1238A and O VI 1032k 

The ubiquity and predominance in network regions of blue- 
ward asymmetries for velocities of order 40-100 km/s in the 
simultaneous C IV 1548 A and Ne VIII 770 A (in second or- 
der) raster scans and the reference spectra in Ne VIII 770 A 
(in first order) are confirmed by several (non- simultaneous) 
SUMER raster scans (taken early during the SOHO mission) 
for three different lines that have high S/N, do not suffer from 
blends and are formed (under equilibrium conditions) at sim- 
ilar temperatures as C IV and Ne VIII. These lines are Si IV 
1402.8 A, N V 1238.8 A and O VI 1031.9 A, formed for 
logT = 4. 8, 5. 3, 5. 5 respectively. 

Figure [9] shows a raster scan with long exposures (60s) of 
the Si IV 1402A line (formed at T ~ 65,000 K) taken at 
quiet Sun at disk center on August 6, 1996. This line does 
not have any blends in the velocity r ange shown (according 
to the SUMER atlas ICurdt et~aT1l2001l) . Similar to the C IV 
and Ne VIII spectra, the R-B analysis reveals that the network 
regions (shown with contours based on core intensity of the 
spectral line) are dominated by red asymmetries for low ve- 
locities, with a sudden shift to blue asymmetries from 50 km/s 
onward. We see significant blue asymmetries even at 80 km/s. 

A very similar picture appears in a disk center raster scan 
(Fig.[[0]) taken about 12 hours later, but now in the N V 1238 
A line (formed at slightly hig her temperatures under equilib- 
rium conditions: ~2xl0 5 K, iMazzotta etaLl 119981) . In the 
network regions, we see a similar shift from redward asym- 
metries for velocities less than 50 km/s, to blueward asym- 
metries for velocities above 50 km/s, up to 80 km/s (the limit 
of the spectral range for this raster). This line has no known 
blends in this velocity range. Again, many locations in the 
network show significant blue asymmetries even at 80 km/s. 

And finally, we consider an R-B analysis of the O VI 
103 1.92 A SUMER observation taken on January 30 1996 
fro m 04:38-04:57UT According to the SUMER spectral at- 
las (ICurdt etal.1 [20011) th e O VI 103 1.9 2 A line (formed in 
equilibrium at ^4xl0 5 K. IMazzotta et al.|[l998l) has one very 
weak blend (< 0.5% of the O VI intensity) in the form of the 
1031. 39A S II line. At this wavelength, any effect of the blend 
in the blue wing will appear at ~154km/s from the line center, 
considerably beyond our analysis window and so we consider 
it to be negligible. Fig. [TT] presents, from left to right, the 
peak line intensity, 1 / e-width of the O VI line for compar- 
ison with a SOHO/MDI line-of- sight magnetogram taken at 
05:21UT, in each case the thick and thin contours delineate 
the brightest network elements and the (approximate) super- 
granular network boundary, through 250 and 100 count peak 
intensity contours derived from panel A. Figure [121 shows the 
results of the R-B analysis of this line, using the steps dis- 
cussed above. As before, we observe that the strong magnetic 
regions that compose the supergranular network exhibit a pro- 
nounced behavior. These regions favor a red-shifted excess of 
material at lower velocities that is replaced by a blue- shifted 
component appearing at about 50km/s with a 5-10% magni- 
tude relative to the peak intensity that persists, gradually fad- 
ing. Therefore, the O VI line has a qualitatively similar be- 
havior to what we have observed above in Si IV, C IV, N V 
and Ne VIII. As an aside, we also notice a very high veloc- 
ity blue component in certain locations (1 12-125km/s - panel 
H). Within the limitations of a slowly rastered image and a 
magnetogram taken 20 minutes later we assess that these lo- 
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Fig. 9. — The August 6 1996 SUMER quiet Sun, disk center, spectrohe- 
liogram in Si IV 1402A (02: 19-04: 18UT - 60s exposure). The top row of 
panels shows the peak line intensity and 1/e line width while the R-B pan- 
els demonstrate the profile asymmetry in steps from 21 to 81 km/s (left to 
right, top to bottom). In each case the quantities are overlaid with inten- 
sity contours at peak intensity levels of 50 and 150 counts to delineate the 
supergranular network boundary (thin solid lines) and the strongest network 
locations (thick solid lines) respectively. The electronic edition of the journal 
has movies showing the panels of this figure. 

cations lie typically (but not exclusively) at the edge of strong 
magnetic elements, are co- spatial with locations of significant 
profile broadening (> 50km/s, see panel B of Fig. [TT]), and 
possibly correspond to "explosive event s" in the periphery of 
the magnetic network (Inne s et aTlll997|) . 

In summary, the raster scans in these three, blend-free, high 
S/N spectral lines show the same blueward asymmetries for 
velocities of order 40-100 km/s in network regions as those we 
found in the C IV and Ne VIII spectra. These results strongly 
suggest that the effect of systemic blends on our results is un- 
likely to cause all observed asymmetries, and that strong up- 
flows of 40-100 km/s seem to occur all over the quiet Sun in 
and around network regions, for lines formed (under equilib- 
rium conditions) at temperatures from the low TR to the upper 
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FIG. 10. — The August 6 1996 SUMER quiet Sun, disk center, spectro- 
heliogram in N V 1238A (14:49-15:19UT - 10s exposure). The top row of 
panels shows the peak line intensity and 1/e line width while the R-B pan- 
els demonstrate the profile asymmetry in steps from 21 to 81 km/s (left to 
right, top to bottom). In each case the quantities are overlaid with inten- 
sity contours at peak intensity levels of 25 and 100 counts to delineate the 
supergranular network boundary (thin solid lines) and the strongest network 
locations (thick solid lines) respectively. The electronic edition of the journal 
has movies showing the panels of this figure. 



A) Peak Intensity |iog l( Counts] B) i/e Lin-? wm-h ;<nV5j ^^Q^tftaggej^jle^QfflSBiL^ 

0.00 0.50 1.C0 1.50 2.00 2.50 3.00 30.00 35.00 40.00 46.00 50.00 -fiS-H>'-16"-W -5 5 10 15 tit S 



S i : ■ I : i '!-! im- H ■ i V ■r:;-i- V- VsjK :-,Oh: :-.-;.iVh-- i : VI ll'.'l.^ ii-'.i h SO' !'-.■■'■/ J I !_■ 2 h "- J - • ■■■■>■■■ •- ''J- 1 - — - 




X (Are Seconds) X (Are Seconds) * lAn; SownOS) 



Fig. 11.— Contextual information for the January 30 1996 1031.93A 
SUMER spectroheliogram of quiet sun at disk center. From left to right 
we show the peak line intensity, the 1/e of the line profile and an 
SOHO/MDI magnetogram taken shortly after the end of the SUMER ob- 
servation (05:21UT). In each case the quantities are overlaid with intensity 
contours at peak intensity levels of 100 and 250 counts to delineate the su- 
pergranular network boundary (thin solid lines) and the strongest network 
locations (thick solid lines) respectively. 
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Fig. 12. — Red-blue (R-B) profile asymmetries (in arbitrary units of inten- 
sity) for the O VI 103 1.93 A emission line. The panels of the figure are over- 
laid with intensity contours at peak intensity levels of 100 and 250 counts 
to delineate the supergranular network boundary (thin solid lines) and the 
strongest network locations (thick solid lines). The electronic edition of the 
journal has a movie showing the panels of this figure along with those at 
intermediate velocities and the components of 
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TR and low corona. These upflows are similar in amplitude 
and nature as those found by D2009 in plage regions, which 
are associated with chromospheric type-II spicules. 

2.4. Time Series ofN V 1238 A 

The raster scans shown in the previous sections clearly 
show a preference for the blueward asymmetries to occur in 
and around network regions. However, all of these scans take 
hours to complete so that temporal characteristics of the up- 
flows cannot be determined. We use a sit-and-stare sequence 
of a quiet Sun region at disk center taken with SUMER on 
25 April 1996 to illustrate the temporal behavior of the up- 
flows. Fig. Q3] shows a one hour long sequence of N V 1238 



A spectra that cover a quiet Sun region, including a network 
region at y=-30 // to y=-20". The intensity of the line is slowly 
changing in the bright network region over the course of the 
one hour long timeseries. The R-B analysis once again shows 
the redward asymmetry for velocities less than 50 km/s which 
shift rapidly to blueward asymmetries for velocities in excess 
of 50 km/s. Blueward asymmetries dominate from 50 to 100 
km/s in the network region. In the internetwork regions very 
little asymmetry is found for these velocities. The right pan- 
els of Fig. [T3l clearly show that the blueward asymmetries are 
highly variable in time, even while they continue to occur in 
and around the network region. This data shows that the up- 
flows of 50-100 km/s are apparently episodic in nature. 

Further analysis of this time series shows that the charac- 
teristic timescale on which these upflow events seem to recur 
is of order 3-15 minutes. This is illustrated in Fig. fT4l which 
shows, as a function of time and space, the episodic nature 
of the R-B asymmetry of N V 1238 A (summed from 58 to 
111 km/s) for the same timeseries. The middle panel shows 
the R-B timeseries at y = — 22" with clear peaks of blueward 
asymmetry occurring, and typical events lasting of order 1-3 
minutes. These events apparently recur on timescales of 3-15 
minutes. This is confirmed by a wavelet analysis (using the 
Morlet wavelet) of the N V 1238A R-B (58-lllkm/s) time- 
series: the resulting wavelet power spectrum (bottom panel) 
shows significant (at 95% significance level) power through- 
out the timeseries at periods of order 3-5 minutes, as well as 
some at 15 min. 

The lifetime of these upflow events agrees surprisingly 
well with those of typ e II spicules (of order 10-100 s 
iDe Pontieu et al.l l2007b|) . This is further supporting evi- 
dence that these quiet Sun events are similar to those seen 
by D2009 in plage regions (which were directly associated 
with type II spicules). It is also very interesting that the 
upflow events in this particular network region occur quasi- 
periodically, or at least seem to occur on typical timescales 
of 3-15 minutes, timescales that are similar to granular (and 
other photospheric) timescales. The quasi-periodicity of the 
R-B asymmetries is further evidence that the periodicities 
observed in coronal imaging and spectral timeseries are not 
necessarily sign s of propagating magneto-acoustic waves i n 
the corona (e.g. Ide Moortel et al.l l2002al: IWang et al.l 12009). 
but may (sometimes) also be caused by quasi-periodic, high- 
speed, faint upflows associated with type II spicules that are 
he ated to coronal temperatures. T his topic is further discussed 
bv lMdntosh & De Pontieul (120091) . 

3. DISCUSSION & INTERPRETATION 

We observe high speed upflows in emission lines that are 
formed, under equilibrium conditions, in the low, middle, and 
upper TR (or low corona). The upflows show a similar spatial 
pattern in all lines, preferentially occur around the network re- 
gions, and have velocities of order 40-100 km/s. The emission 
of the upflowing plasma is faint at ~ 5% of the peak intensity 
of the line. It is possible that th ese upflows are visib le as 
propagating disturbances in EUV ([Sc hriiver et al. 1999) and 
soft X-Ray coronal images (Sa kao et al.ll2008|) : the velocity 
of the measured weak episodic jets are in the same range as 
the apparent velocities of the disturbances in imaging data. 
A connection to the chromospheric activity and/or spectro- 
scopi c measurements is studied by (Mcintos h & De Pontieul 
2009) for active regions. 

What causes these upflows? Many authors have dis- 
cussed blueward excursions of TR lines in the past 
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Fig. 13.— The April 25 1996 SUMER quiet Sun, disk center, timeseries in N V 1238A (22:01-23:01UT - 7.5s exposure). The panels, from left to right, show 
the peak line intensity and the R-B profile asymmetry in steps from 15 to 105 km/s. The panels are overlaid with intensity contours at peak intensity levels of 
75 and 150 counts to delineate the supergranular network boundary (thin solid lines) and the strongest network locations (thick solid lines) respectively. The 
electronic edition of the journal has movies showing the panels of this figure. 
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Fig. 14.— Wavelet analysis of the N V 1238A R-B (58-lllkm/s) time- 
series (top panel, see e.g, Fig. |T3). The R-B timeseries at y = —22" are 
selected for wavelet analysis (center panel) and the resulting wavelet power 
spectrum is shown (bottom panel). The bottom panel shows the cross-hatched 
cone of influence and contours of the 90% significance level in the timeseries. 



(e.g. iDere. Bartoe & Brueckner 1989]), with most of those 
observations being linked to explosive events: vio- 
lent, often bi-directional jets caused by reconnection 
(Chae, Schiihl e & Lema ire 1998). What we report on here is 
not an occasional disturbance, like an explosive event, but a 
much more ubiquitous and fainter phenomenon: at any sin- 
gle moment, most of the network regions show these high 
speed upflows. We believe that these upflows are related to 
the blueward asymmetry and line broadening that has been 
reported b efore for TR lines in the quiet Sun network (e.g. 
lPeterll2Q()Qh . The latter suggested, without the benefit of the 
R-B measure, that hot coronal funnels cause this faint compo- 
nent, with the dominant core of the line emitted in small-scale, 
low-lying coronal loops. Our interpretation is very different. 
The quiet Sun upflows are very similar in location, velocity 
(50-120 km/s), intensity (5% of line core) and temporal be- 
havior to the faint, coronal upflows observed with Hinode/EIS 



above active region plage (lHara et al.ll2008l D2009). The lat- 
ter have been directly linked to the heating to coronal temper- 
atures of chromospheric type II spicules, which dominate the 
upper chromosphere in plage and quiet Sun network and de- 
velop velocities of order 50-150 km/s. D2009 estimated that 
the mass flux injected to the corona by these jets is significant 
for the mass balance of active region coronal loops. We sug- 
gest that the upflows we observe here in Si IV, C IV, N V, 
O VI, and Ne VIII are the quiet Sun equivalent of this process. 

In Fig. [15] we provide a cartoon of the TR structure consis- 
tent with our observations and the analysis of D2009. What 
we observe as the transition region is a conveyor belt of sorts, 
one that includes frequent bursts of rapidly upflowing ionized 
hot plasma (in the form of thin, dynamic type II spicules), 
and recycles the material that is slowly cooling out of the 
corona, losing energy through radiation and producing the 
well known prevalent TR red- shift, as we ll as a host of other 
TR phenomenology discussed above (e.g. , iGebbie et al.ll 198 lb 
de Wijn & De Pontieu 2006; Judge 2008). This figure is very 
different from Fig. 8 of Peter ( 2000). We interpret the super- 
granular magnetic network as a truly dynamic environment of 
very high speed ionized plasma upflow (in the form of type II 
spicules) with a mix of instantaneous and gradual mass return 
from the corona instead of a component due to "hot coro- 
nal funnels" rooted in su pergranular boundaries. While low- 
lying, cooler loops (e.g.. lDowdy. Rab in & Moore 1986) may 
well play some role 8 in the TR emission (lPeterl2 000) our anal- 
ysis has demonstrated that the role of spicules is fundamental 
in understanding the observed TR signature. It is possible 
that the small network core blueshifts observed in Ne VIII 
may be the result of a more gentle evapo rative mass load- 
ing due to coronal energy d eposition (e.g., [Klimchuk 2006; 
Tomczyk & Mcintosh 2009) in the upper TR. Regardless of 
the interpretation, our results suggest that approximating the 
observed spectral line profiles with (even multiple) Gaussians 
is profoundly limiting, and that other ways to consistently 
model the observed asymmetric p lasma velocity distributions 
are essential (e.g., Scudder 1992a|). 

Our c onceptual view of t he TR is related to the model put 
forth by Pneu man & Koppl (Q978) who suggested that a small 
fraction of spicules are heated to coronal temperatures and 



8 Judge & Centeno ( 2008) recently questioned this view of transition re- 
gion structure. 
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Fig. 15. — Conceptual view of TR emission in network and inter-network. 

provide the corona with hot plasma. To constrain the spicular 
mass flux into the quiet Sun corona we follow the approach of 
D2009, appealing to the conservation of mass and observed 
emission measures in the system and noting that the latter are 
similar to those observed in active regions - upflow compo- 
nents with ~ 5% of the core intensity. D2009 found that 
the density of the coronal component of the spicules could 
be expressed as pi « p c / (3e^) where p c the coronal den- 
sity, and €t is the fraction of total coronal spicular density 
at a temperature of T. If we assume that only a fraction 
ct = 0.25 — 0.5 of the spicular coronal mass flux is observed 
in Ne VIII, then the only difference with the analysis of D2009 
is that the coronal density in quiet Sun is less than that of ac- 
tive regions (say, 3xl0 8 instead of 10 9 cm -3 ). As a result 
the coronal densities p c of hot spicular matter are of order 
3 — 15xl0 8 cm -3 . This implies that the heated spicules can 
fill the quiet Sun corona even if only 0.3-1.5% of the chromo- 
spheric spicule mass reaches coronal temperatures. D2009 
also estimate that the number of spicules per arcsec 2 oc- 
curring within the network over the lifetime of a spicule is 
N n « 2e|725 2 /(5 2 , where Si « 250km the diameter of a 
typical spicule. Since the assumptions about the coronal scale 
height, expansion from network to corona, velocity and diam- 
eter of spicules are of the same order for active region and 
quiet Sun, we estimate a density of spicules of order 1-4 per 
arcsec 2 within a network region. Is this number compatible 
with Hinode/SOT observations of type II spicules at the limb? 
IRouppe van der Voort et al.l (|2009|) report that they see at least 
between 1 and 3 of these s picules per linear arcsec at th e limb 
(at a height of 5,000 km). IRouppe van der Voort etafl (120091) 
indicate that when observing at 5,000 km height above the 
limb, line-of-sight superposition implies that one can see fea- 
tures emanating from network regions over a line-of-sight that 
is of order 200 arcseconds long. Such a line-of-sight in prin- 
ciple could cross network at the edge of of ~ 4 supergranules. 
However, the network around supergranular cells has a rela- 
ti yely low filling factor of orde r 0.25, as shown in, e.g., Fig. 12 
of Roupp e van der Voort et al.l (|2009). This means that our es- 
timate of 1-4 spicules per arcsec 2 within a network region fits 
in well with the observed number of spicules (1-3 per linear 
arcsec) that is observed with Hinode at the limb (since the 
number of supergranule cells multiplier is cancelled by the 
low network filling factor). We conclude that the observed 
number of spicules at the limb is fully compatible with a sce- 
nario in which the heated component of type II spicules ob- 
served in TR emission can play a significant role in the mass 
balance of the quiet corona. 



The impact of these spicules on the coronal energy balance 
can be estimated by calculating the flux associated with the 
internal energy of the mass flux calculated in the above. We 
find that the energy flux is given by Nip c c v Tvi, with Ni the 
number of spicules that occur at any location over the spicule 
lifetime, p c the mass density of the spicule component that is 
heated to coronal temperatures, c v the specific heat at constant 
volume, Vi the upflow velocities in the spicules, and T the 
temperature the coronal spicular component reaches. D2009 
estimated that Ni « 2ef,. If we assume that all the plasma 
in these spicules is heated to 1 MK (of order the formation 
temperature of Ne VIII, this implies er ~ 1), and use p c = 
3 x 10 8 cm— 3, Vi ~ 100 km s -1 , we find an energy flux 
of order 7 x 10 5 erg cm -2 s -1 . This is of the same order 
of magnitude as the estimated conduc tive and radi ative losses 
(3 x 10 5 erg cm -2 s _1 for quiet Sun. lPriestl[i 9821) . 

Several issues are unresolved in the scenario presented here. 
The first is the role that spicules play in the mass balance of 
the lower TR. Observations from Hinode/SOT indicate tha t 
many are heated above 20,000K (iDe Pontieu et al.l I2007bl) . 
but the peak temp erature that the bu l k of th e spicular mass 
reaches is unclear. iPneuman & Kopp (119781) suggest that the 
fraction of spicular mass reaching a certain TR temperature 
T is a steeply declining function of T, consistent with the 
rapid decrease of differential emission measure from 10,000 
to 100,000K. Simultaneous observations with high spectral, 
temporal and spatial resolution of the chromosphere and TR 
at high signal-to-noise (>20) are necessary to resolve this is- 
sue. With the recent selection of the Interface Region Imaging 
Spectrometer (IRIS) such data should be forthcoming. 

A second issue relates to the formation and heating mech- 
anism of type II spicules themselves. Insight into the mag- 
netic geometry required for the energy release and how the 
plasma itself is heated and accelerated along the structure of 
the spicule may help explain why we see a velocity differ- 
ence in the upward velocity components between the forma- 
tion conditions of Si IV, C IV, N V, O VI and Ne VIII as 
well as those between the quiet and active atmosphere (with 
upflows stronger in the latter - up to 150 km/s). Certainly, the 
temporal behavior points to a driver that leads to episodic up- 
flows on ch aracteristic timescales of granulati on, something 
also seen by Roup pevan der Voort et al.l (|2009) who report on 
the discovery of the disk-counterpart of type II spicules. 

Finally, based on the observational evidence and interpreta- 
tion presented in this paper, efforts to better understand type 
II spicules may lead to a clearer picture of the impact of non- 
equilibrium ionization on the upper portion of the solar atmo- 
sphere. 
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